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Abstract 

We analyze the kinetic behavior of localized excitations - solitons, breathers and 
phonons - in Sine-Gordon model. Collision integrals for all type of localized excitation 
collision processes are constructed, and the kinetic equations are derived. We prove 
that the entropy production in the system of localized excitations takes place only in 
the case of inhomogeneous distribution of these excitations in real and phase spaces. 
We derive transport equations for soliton and breather densities, temperatures and 
mean velocities i.e. show that collisions of localized excitations lead to creation of 
diffusion, thermoconductivity and intrinsic friction processes. The diffusion coefficients 
for solitons and breathers, describing the diffusion processes in real and phase spaces, 
are calculated. It is shown that diffusion processes in real space are much faster than 
the diffusion processes in phase space. 
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1. INTRODUCTION 



The problem of kinetic properties of excitations in integrable models belongs to a 
class of most nontrivial problems of physical kinetics. First, enormously long relax- 
ation of nonlinear excitation has been found in numerical experiments by Fermi Pasta 
and Ulam [1]. In Zabusky and Kruskal numerical experiment [2] unexpected behavior 
of localized nonlinear excitations has been discovered, namely their interaction with- 
out changing their forms and velocities. Zabusky and Kruskal named them solitons. 
Shortly after, the analytical method of solving nonlinear differential equation with par- 
tial derivatives - inverse scattering method - was found in the framework of Korteweg 
- de Vries equation for 1+1 dimensions [3]. Other physically reasonable continuous 
models, treatable by the inverse scattering method are the Nonlinear Schrodinger (NS) 
equation, the Sine-Gordon (SG) equation and the Landau-Lifshitz (LL) equation [4,5]. 
Common excitations in integrable models are: the one parameter localized wave i.e. 
the soliton and the nonlinear periodic wave for KdV model; and the two parameter 
localized wave i.e. the breather for SG and LL models; and the envelope soliton for 
NS equation. All these excitations interact without changing their forms and their 
velocities and therefore their energies; the only result of their interaction are shifts 
of their coordinates and for breathers and NS solitons change of their phases as well. 
It is worthwhile to emphasize that many particle effects are absent in the following 
meaning: the total shift in collisions involving several excitations is equal to the sum 
of shifts in each collision (see[4]). 

Soon after analytically solving the KdV equation [3] , Zakharov considered the pos- 
sibility of the kinetic decryption of solitons [6]. In [6], kinetic equation for solitons 
was written but only the effect of renormalization of soliton velocity due to soliton - 
soliton collisions has been considered. Another approach for the investigation of ki- 
netic properties of kink type solitons was proposed in [7-9]. In [7] numerically and 
in [8,9] analytically the diffusion of kinks interacting with phonons in the 4 model 
was considered. This model is not exactly integrable, but the potential for kink - 
phonons interaction has reflectionless form and the result of the interaction is the same 
as for soliton - soliton collision in integrable models. Taking into account such kind 
of interaction in [8,9], kink diffusion coefficient was calculated. A bit later the same 
calculations have been done for SG model [10]. The renormalization of soliton velocity 
and the diffusion coefficient of the solitons due to soliton - soliton and soliton - magnon 
collisions have been considered for SG model in [11]. 

To explain the effect of velocity renormalization it is useful to note that the depen- 
dence of soliton coordinate x from time can be written in the following form: 



x(t) = x + v\t + / dt' I dx 2 dv 2 \v\ — v 2 \Ax s (v 1 ,v 2 ) 6 (t — t') f(x 2 ,v 2 ,t'), (1.1) 



when shifts of soliton coordinates Ax s (vi,v 2 ) are taken into consideration. In formula 
(1.1) f(x,v,t) is the distribution function of solitons or phonons (magnons), and the 
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sample soliton collides with soliton or phonon (magnon) wave packet (#2, 1*2) at the 
moment t', 

e(t) = i, t> o, o(t) = o, t< 0. (1.2) 

^From formula (1.1) one can obtain the following expression for the average velocity 
of the sample soliton: 

< vi >= vi + J dx 2 dv 2 \vi — v 2 \Ax s (v!, v 2 )f(x2,v 2 ,t) (1.3) 

The diffusion coefficient of soliton is give by: 

D _ <[x(t)-<v>t-x(0)} 2 > ^ 

(see [8-10]). 

So the shift of soliton coordinate Ax leads to its diffusion. Do another kinetic 
coefficients for solitons exist in the integrable models? This question is closely related 
with the problem of entropy production. ^From the fact of conserving of the solitons 
velocities follows the absence of the entropy production in momentum space only. But 
in the coordinate space the interaction of solitons is very strong. Therefore the question: 
" Does entropy increase in integrable systems?"- is not a nonsense. The answer on these 
questions have been done in [12] on the example of kink-type solitons in SG model. In 
[12] kinetic equation for the gas of kinks has been formulated, the entropy production 
has been proved and the existence of the coefficients of intrinsic friction and thermal 
conductivity have been shown. Obviously that such kinetic coefficients as mobility 
coefficient appear in the presence of perturbations, which destroy the integrability of 
the model ( see [13 ] and review [14] ). 

In this paper the kinetic behavior of solitons, breathers and phonons in the frame- 
work of the SG model is considered. We construct collision integrals for all possible 
collisions and formulate the system of Boltzmann type kinetic equations for solitons, 
breathers and phonons. Based on the system of kinetic equations thus obtained, we 
prove that the entropy production takes place as a result of randomization of the 
distribution of excitations in coordinate and phase space. Thus we are able to derive 
transport equations for solitons and breathers and to calculate self-diffusion coefficients 
for soliton-soliton and soliton-breather and breather-breather collisions as an example. 

2. ELEMENTARY EXCITATIONS IN SINE -GORDON EQUATION 



Sine-Gordon equation in dimensionless variables can be written as: 



2 

m 



<Ptt ~ ( Pxx + ^- sin f3(f = 0. (2.1) 
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Here ip is the order parameter, t is the dimensionless time, x is the dimensionless space 
coordinate, c=l is the characteristic velocity, m is the dimensionless mass and (5 is the 
parameter of nonlinearity. 

Eq.(2.1) follows from Hamilton's equations (see, for example, [5]): 



ip t = {H, (fi},7r t = {H, %} 



(2.2) 



with Hamiltonian: 



2 



m 



7T 2 + ^ + 2l-l (l-COS/fy>) 



dx 



(2.3) 



In eqs. (2.2) -(2.3) ip and n- are canonically conjugated coordinate and momentum: 

{n(x),<p(y)} = 6(x-y), (2.4) 
where the Poisson brackets are defined by usual way: 

' 5A 5B 5A 5B 



{AB} 



dx. 



k 5tt(x) Scp(x) Sip(x) 87r(x) t 
It is easy to see that the functionals of total momentum, P, of the system: 



(2.5) 



P = -J n(x,t) 
and total "angular" momentum K: 



dip (x, t) 



dx 



dx 



K 



J xh [tt, ip] 



dx 



(2.6) 



(2.7) 



commute with Hamiltonian (2.3). 

It is well known that there are two types of localized excitations (LE) in SG model: 
solitons and breathers. Soliton is a one parameter solution, breather is a two parameter 
solution, the breather can be interpreted as bound states of two solitons (see [4,5]). 
The one parameter solution of the SG model, often named as kink, can be written as: 



P 

Breather solution has the following form: 



ip s = -A^arctg exp (S s 1 (x - v s t - x 0s ) 



f u 2 \ sin [u(v b )t - k(v b )x - ip ob ] 

ip b = -arctg — r/ yr-, 

V^i / ch lo b (x — v b t — x ob 1 
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(2.9) 
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where: 

m(v) = m/yl — v 2 , u(v) = m(v)ui, k(v) = vm(v)ui, 

8? = m(v a ), 5, 1 = u 2 m(v b ) . (2.10) 

Here e={+l,-l,0} is the soliton topological charge, +1 corresponds to soliton, -1 to 
antisoliton, to breather, v s and Vb are the soliton and breather velocities, 5 S and 5b 
are the soliton and breather sizes, xo s and xob are the initial coordinates of soliton 
and breather, (/?o6 is the initial breather phase, and u> 2 are the breather parameters 
satisfying the following condition (see [5]): 

u\ + ul = \. (2.11) 

When u 2 — * 0, the breather solution reduces to the following ones: 

hm .(pb(x,t) = (p p (x,t) = 4— — ? — rr, 2.12 

and gives a plane wave solution i.e. a phonon. 

The phonon frequency is related with wave vector k by formula: 

u 2 ph (k) = m 2 + k\ (2.13) 

The breather velocity, when oj 2 — * 0, reduces to the phonon group velocity: 

v ph = (k/ujp h ) = (dujp h /dk) , (2.14) 

and 4u 2 becomes the amplitude of phonon oscillations. 

^From formulas (2.8), (2.9) and (2.3) it is easy to find the soliton and the breather 
energy: 

E s = M a /y/l - vl E b = M b /y/l-vl (2-15) 
The soliton and the breather rest masses are given by: 

M s = 8m, M b = 16muj 2 . (2.16) 

When uj 2 — > 1, the breather reduces to soliton-antisoliton bound state with a mass 
16m. The difference between the breather energy on and the sum of energies of isolated 
soliton and antisoliton gives the binding energy of the breather: 

16?7l 

A£ = -- 7 _(l-c 2 ) (2.17) 
If the system consists of n\ solitons and n 2 breather we have: 

ni ni ni «2 ni ni 

E = Y. E s + Y.E b , P = J2P s + Y / P b , K = J2Ks + Y / K b . (2.18) 

s=l 6=1 s=l 6=1 s=l 6=1 
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Energy and momentum are related by following formula: 



E 2 = M 2 + P 2 . (2.19) 

Formulae (2.18) allow us to consider solitons and breathers as elementary excitations 
of field (p. 

3. COLLISIONS OF LOCALIZED EXCITATIONS 

Due to the ID of the problem the collision of solitons and breather take place without 
dependence from values xo s and xob , s — 1, ni, b — 1, n 2 because their velocities do 
not change. This peculiar property follows from energy and momentum conservation 
laws in Sine-Gordon equation as well as other exactly integrable models (see [4,5]). 
Those conservation laws have the form: 

vi = v[; v 2 = v' 2 . (3.1) 

Here and later the values after collision will be devoted with a prime. 

To analyze collision processes it is necessary to take into account the the " angular" 
momentum conservation law, which can be written as: 

x x E x + x 2 E 2 = x[E[ + x' 2 E' 2 . (3.2) 

The result of collision is coordinate shifts (changes) and for breathers phase shifts as 
well. Furthermore only pair collisions exist and there are no many particles effects [4]. 
Therefore formulae (3.1), (3. 2) provide the general framework for studying the effects 
of collisions. In this section conservation laws and coordinate Xqi and phase changes 
(fioi for all type of collisions will be written down explicitly from general formulas [4,5]. 
For simplicity the index "0" will be omitted. 

Soliton-soliton collisions. Conservation laws for two soliton collision have the fol- 
lowing form: 

v' X8 = vis, v' 2s = v 2s , (x' ls - x ls )E ls + (x 2s - x 2s )E 2s = (3.3) 
Coordinates shifts are given by: 

4 Sis 
Axis = x'is - x u = -E^sign(v ) In \Z SS \ = -^-sign(v ) In \Z SS \ 

tils I 



Ax 2s = x 2s -x 2s = -— sign(v ) In \Z SS \ = — ^-sign(v )\n\Z ss \, (3.4) 



where 



1 + yjl-vt {y ls - y 2s ) 

Z ss = i = , = Ti V- 3 - 5 
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Here vq is the velocity of relative motion of solitons which in the general case can be 
written as: 

^o = 7T^ — ~~~~T? i,k = s,b,ph. (3.6) 

(1 - VuV2k) 

Soliton - breather collision. In this case conservation laws can be written as: 

vu = v' u , v 2b = v' 2b , u 22 = uj' 22 , 

(x' ls - x u )E u + (x' 2b - x 2b )E 2b = 0. (3.7) 
Coordinate shifts Ax are described by the following formulae: 

g 

Ax ls = x' ls - x ls = —sign(v ) In \Z sb \ = 8 ls sign(v ) In \Z sb \ 
E\ s 

Ax 2b = x' 2b - x 2b = ~-^-sign(v ) In \Z SS \ = -^-sign(v ) In \Z sb \, (3.8) 

&2 b £ 

where: 



1 + ^2, b 2\jl-vl ( Vls - V 2b ) 

Vo = j- r. (3.9) 



1-^2,62^1-^0 0--VnV2h) 

The breather phase shift (Atp) b is determined by the formula: 



tan(A</?) 6 = —sign(vo)- ° ^ . (3.10) 

Some comments regarding formulae (3.8) and (3.10) are in order. For soliton coordinate 
shift (formula (3.8)) due to soliton - breather collision one has a factor of "8" instead 
of "4" as in formula (3.4). This difference can be explained as follows. If uj 2tb2 — > 1, one 
can considered a breather as two solitons with opposite topological charge. Since the 
coordinate shift does not depend on topological charge, the soliton - breather collision 
in this case can be considered as a soliton - two solitons collision, and thus the shift is 
doubled. Let us mention that in the case cu 2 ,b2 1, the quantities Z ss and Z sb are the 
same. 

Formula (3.10) implies that the phase shift is in general large: {A(p) b ~ tt. It 
is small one only in the ultrarelativistic case \v s — 1| << 1. In this case (A(p) b i pa 



-sign(yo)\jl - i>o^i,&i -If v s « 1, then (A(p) bl pa tt/2. 

Soliton-phonons collision. This process is characterized by the following conserva- 
tion laws: 

V'u = V ls , V 2ph = V 2ph , J ph = LJ ph , 

(x' ls - x ls )E ls + (x 2ph - x 2ph )E 2ph = 0. (3.11) 
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The expressions for x and cp shifts are: 



mJl — v 2 s 

Ax 2ph = x' 2ph - x 2ph = swn{v Q ) uA ^ h _ kVs) (3.12) 



Axi s = x[ s - xis = -^r^Az2 P /i (3.13) 

Jl-vl 

tan(Ay9 2p/l ) = -sign(v )- . (3.14) 

Phonon energy E 2ph and frequency cu p h(k) are defined by the formulae: 

E 2ph = 16u 2 u ph (k), u 2 ph (k) = m 2 + k 2 . (3.15) 
Breather - breather collisions. The corresponding conservation law have the form: 

v'u, = V U , V 2b = V 2b , UJ 2 bl = UJ 2 ,bl: ^2,62 = U 2,b2 

(x' lb - x lb )E lb + (x' 2b - x 2b )E 2b = 0. (3.16) 
The expressions for the x b and ip shifts are: 

Ax lb = -^-sign(v ) In \Z bb Z' bb \ = ^-sign(v ) In \Z bb Z' bb \ 

Ax 2b = -Jt-sign(v ) In \Z bb Z' bb \ = -^-sign(v ) \n\Z bb Z' bb \ (3.17) 

&2b * 

2vq v/l — Vq sin tpi cos ip 2 

tan(A^ lfe ) = sign(v ) 2 v — ^- — — -. (3.18) 

Vq + (1 — Vq) (cos 2 ipi — COS 2 



2f \Jl — v 2 . sin ^ 2 cos ^1 

Vq + (1 — t>o)(cOS 2 ^1 — COS 2 ^2) 



tan(Av? 2 fe) = -sign{v ) -^ — 2 T77Z3. ~ ( 3 - 19 ) 



Here: 



1-^1- Ugc08(^l+ -02) , l + ^/l -ugcOS^l --02) 

^66 — 1 ) ^ bb — 1 • {6.2,0) 

l-y 1 -uocos^i -^2) 1 + V 1 -^o cos (^i +^2) 

The angles ip\ and -0 2 are related to the parameters cj 2 ,6i and u) 2 , b2 by the formulae: 

sin^i = uj 2yb i, sin ^ 2 = ^2,62 ( 3 -21) 

Breather - phonon collision. Since phonons are the limiting case of breather when 
uj 2 — > 0, formulae for breather - phonons collision can be obtained from ones in the 
previous paragraph by passing to the limit a; 2 62 — > 0. It means that the second breather 



8 



reduces to phonon wave packet and the notion of group velocity and coordinate of wave 
packet automatically appears. The breather velocity v b reduces to the group velocity 
of the wave packet and the coordinate x b reduces to the coordinate of the center of the 
wave packet. 

The corresponding conservation law has the form: 

V' lb = Vlb, V' 2ph = V 2ph , W' 2j bl = W 2,W, ^2,ph2 = W 2,ph2 

(x' lb - x lb )E lb + (x' 2ph - x 2ph )E 2ph = 0, (3.22) 

where: 



2muj 2)b iJl — v b E 2p h 

Ax 2ph = -sign(v ) ; : — ^ T 2 . 2 1 , Ax xb = — —Ax 2ph (3.23) 

u P h[u P h - kv b ][l - (1 - v§)uf tbl ] E lb 

2v Jl — VqUJ 2 5! 

tan(Av? 16 ) = -sign(v )— 2 f— ' (3.24) 

ZV — 1 + {I — V )LU lbl 



a . i v 2v \/l-v$u 2>ph2 u 1M 

A^ 2ph = -sign(v ) . (3.25) 

zv — 1 -t- {1 — v )u l bl 

Phonon-phonon collisions. As it is well known, in linear theory phonon-phonon 
collisions are absent. Here phonon - phonon interaction is the result of nonlinearity 
of the Sine-Gordon equation. Corresponding relations can be obtained from ones in 
previous paragraph by passing the limit u 2tb i — > 0. By considering cu 2>b i as a small 
quantity, cu 2>b i << 1, one has: 

v i P h = v iph, v' 2ph = v 2ph , Axi ph Ei ph + Ax 2ph E 2ph = 0. (3.26) 

The shift of phonons coordinate is: 



Axi ph = sign{v ) — — 2 Ax 2ph = -— — Ax lph , (3.27) 

£>l,ph v &2ph 

where: 

Eiph — ^G^iph^iph, E 2ph = 16uj 22p hUJ 2 ph- (3.28) 
The change of phonon phase is: 



Aifip h = -2sign(v )u 2 , phl , A<p 2ph = -2sign(v )uj 2tPh2 . (3.29) 

Vq Vq 

The examination of phonons as a limiting case of breathers gives the possibility to 
consider changes of phonon phase and coordinate simultaneously. 
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Let us analyze phase shifts of breathers and phonons. It is easy to see from formulae 
(3.10)-(3.12) and (3.14) that breather and phonon shifts are small for soliton-breather 
and soliton-phonon collisions in the general case. For breather - phonon collisions 
breather and phonon shifts are essentially different (see (3. 25), (3. 29)). Breather phase 
shift are proportional to to 2p h an d therefore small. Phonon phase shift (Atp) p f l ~ 1. For 
phonon - phonon collisions phase shifts are small (see (3.33)). This agrees with the 
standard linear theory of phonons. In accordance with the usual theory of phonons 
shifts of phonon wave packets do not take place. 

In conclusion, let us emphasize that soliton topological charge is conserved for each 
type of collisions. The processes leading to soliton-antisoliton bound state formation 
and the reverse processes do not take place because of energy conservation. Thus in 
the framework of the exactly integrable Sine-Gordon model with Hamiltonian (2.3) 
the numbers of solitons, antisolitons, breathers and phonons are defined by initial 
conditions at each point of (x, t). Interactions destroying the integrability of the system 
lead to possibility of energy exchange between the LE and to the formation and decay 
of breathers. In the present paper only the processes with conservation of energy and 
topological charge have been considered. The processes with energy exchange will be 
analyzed in a future publication. 

4. PROBABILITY OF SCATTERING PROCESSES 



To construct collision integrals for processes considered in the present section let us 
analyze the probability of the corresponding scattering process. The simplest process 
is soliton - soliton scattering. By definition the initial state (i) of a soliton pair is 
described by its velocities and center of gravity coordinates: 

i = (x 1 ,v 1 ,x 2 ,v 2 ) = (1,2) . (4.1) 
The final state (f) is defined by: 

/=K,^,4^ 2 ) = (1',2'). (4.2) 

Current density of solitons per unit density is given by: 

Ji = \vo\- (4.3) 

Taking into account formula (4.3) and the fact that two solitons collide in any case, 
the transition probability per unit time from state (i) to state (/) can be presented as: 

W w = W(l', 2'|1, 2) = \v \5(v[ - Vl )8(v' 2 - v 2 )5&x' 1 - Xl )+ 

hi 2 

+(x 2 - x 2 ))5(x[ -xi- Ax(ui, v 2 )). (4.4) 
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In this formula the first two delta functions describe energy conservation laws, the third 
one - " angular" momentum conservation law, and last delta function describes coordi- 
nate shift. The coordinate shift of second soliton follows from the "angular" momentum 
conservation law. The coordinate shift A x\ is defined by expression (3. 4), (3. 5). 
It is convenient to rewrite formula (4.4) in the form: 

W(l',2'|l,2) = 

= R ss (v' 1 ,v 2 \v 1 ,v 2 )5(x' 1 -xi- Ax 1 (v 1 ,v 2 ))S(x' 2 - x 2 - Ax 2 (v 1 ,v 2 )). (4.5) 

where: 

R S a(v' 1 ,v' 2 \v 1 ,v 2 ) = \v \5(v[ - vi)S(v' 2 - v 2 ). (4.6) 

Let us note that: 

RssW^v'^v^v^ = Rs S (vi,v 2 \v[,v 2 ) = Rss^v'^v^vx). (4.7) 
For the probability of the inverse process the following expression can be written: 

W(l,2|l',2') = 

= R ss (v' 1 ,v 2 \v 1 ,v 2 )5(x 1 - x[ + Ax 1 {v 1 ,v 2 ))8(x 2 - x' 2 + Ax 2 {vi,v 2 )). (4.8) 

Let us emphasize that the probability defined by formula (4.4) does not satisfy the 
detailed balance principle in the standard form: 

W{1\ 2'|1, 2) = W(V, 2*j 1'*, 2'*). (4.9) 

In (4.9) the following notation has been used: (1*) = (xi, — i>i) . 
For further analysis it is convenient to present the expression (4.7) for the proba- 
bility W as a sum of two parts W r and W m : 

W ss (l, 2|1', 2') = W; a (l, 2|1', 2') + W™(1, 2|1', 2'), (4.10) 

where: ^ 

W L = -R ss [5(x ls - x' u - Axi s )5{x 2s - x 2s - Ax 2s )- 

-5(x ls - x' u + Ax ls )5(x 2s - x' 2s + Ax 2s )] (4.11) 

w Ts = \Rss[K x is ~ x' ls - Ax ls )5{x 2s - x' 2s - Ax 2s ) + 

+5(x ls - x' u + Ax ls )5(x 2s - x' 2s + Ax 2s )]. (4.12) 
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The probability WJ S , as it will be shown later, describes the effect of renormalization 
of soliton velocity, and the probability W™ describes the homogenization of solitons 
phenomena i.e. the effect of mixing of states. The quantity W™ is related with entropy 
production in the soliton gas and the kinetic coefficients (see sections 5,6). 
It is easy to convince oneself that: 

W^(1 , ,2'|1,2) = W^(1,2|1',2'). (4.13) 

This equality means that for the dissipative part of the transition probability the usual 
detailed balance principle is valid. 

Let us consider now the general case of LE i and k, (i,k = s,b,ph) for soliton, 
breather and phonon correspondingly. The probability of such type of collisions per 
unit time W ik can be defined by following formula: 

W ifc = Wi fc (l , ,2 , |l,2) = 

= R lk 5(X[ -X % - d t (l, 2))5{X' k -X k - 4(1, 2)), (4.14) 

where: 

Rik = Rik{V!, V k '\V u V k ) = \v \5(Vi - Vf)6(V k - VI). (4.15) 

In formulae (4. 14), (4. 15) the following notations have been used. Numbers 1 and 2 
mean the set of variables defining of the colliding LE states, dj, is the shift of LE 
coordinate or phase. In the case of solitons 1 = (x s ,v s ), for breather or phonon 
1 = (xi,ifi,Vi,u)2i), i = b,ph. For simplicity two component coordinate Xi and two 
component velocity Vi have been introduced in the following way: 

X s = x s ,V s =v s , Xi = (xi,ipi),Vi = (vi,u 2 i),i = b,ph. (4.16) 

In (4.15), the delta function, 8 {Vi — V/), describes conservation law of energy and 
momentum of each LE; in the expression (4.14), the delta function 5 (Xl — Xi — d^l, 2)) 
defines the coordinate and phase for breather and phonon shifts after collision. 
^From the definition of Ri k it follows that : 

Ri k = R(V;, VlW, V k ) = R(V, V k \Vj, Vi). (4.17) 

Later only the case of weak inhomogeneous gas of LE will be considered, i.e. when the 
characteristic length / of the distribution function is much larger than the shift AXf 

\h\ » \AXi\ (4.18) 

The mean free path of LE is: 

/ ~ -, (4.20) 

n 

because the cross section of scattering is equal to 1. 
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The coordinate shifts of LE are of the order of Si (see section 3), where 

6i(v,u 2 )~—±— (4.22) 

Therefore the condition (4.18) rewritten in the form: 
1 1 

— >> Si or m(v )u2 >> j ~ rii i — s,b,ph, (4-21) 

Tli I 

means that average distance between LE much larger their sizes. In other words the 
condition of small inhomogeneous (4.18) coincides with the condition of small density 
of LE gas. 

Let us present Wik as a sum of two parts W[ k and by the same way as (4.10)- 
(4.12). Assuming that (4.18) is valid, we can expand the delta functions in the expres- 
sion (4.14) for Wifc(l,2|l,2) in powers series AAj. Then the expressions for W[ k and 
W% become: 

W^(l',2'|l,2) = R ik (l'2' : \l : 2)[di^- + d k ^-}S(Xi-X^S(X k -X' k ) (4.22) 

dXi dX k 

W£(l',2'|l,2) = 

= R ik {l'2', |1, 2) [2 + + d\^ + 2did k -^^-]S(Xi - X')S(X k - X' k ). (4.23) 



5. COLLISION INTEGRALS 



In the present section the expressions for the collision integrals for all type of LE in 
SG equation will be derived. It is easier to do this starting from the formulae for 
probabilities of collisions. Let us introduce following notations: 

f = f(X,V,t),B = B(X,V,t),N = N(X, V, t). (5.1) 

for the distribution functions of solitons, breathers and phonons as / ; B and iV corre- 
spondingly. The probability to find the LE in the state (1, 1 + dl) can be defined the 
usual way: 



dWi = Fi(l)dl, (5.2) 

where 

% = s,b,ph and F s = f,Fb = B,F p h = N. The total result of collisions of a sample 
LE with the other LE is the sum of each collision. This is due to the special type of 
interaction in exactly integrable models. Therefore the general collision integral can 
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be presented as a sum of partial collision integrals. For example in the case of solitons 
the collision integral can be written as: 

C s = C ss {f, /} + C sb {f, B} + C sph {f, N}, (5.3) 

where £ ss , £ s b, £ sp h are soliton-soliton, soliton-breather and soliton-phonon collision 
integrals respectively, which have the following form: 

C„{f,f} = J dl , d2 / d2{W;.(l,2|^2 / )/(lO/(20-W M (l / ,2 / |l,2)/(l)/(2)} (5.4) 

Csb{f,B} = J dl'd2'd2{W sb (l,2\l',2')f(l f )B(2')- 

-W sb (l',2'\l,2)f(l)B(2)} (5.5) 

£ sph {f,N} = J dl'd2'd2{W sph (l,2\l',2')f(l')N(2')- 

-W ss (l',2'\l,2)f(l)N(2)}. (5.6) 

The first term in formula (5.4) and in formulae (5.5) and (5.6) describes solitons "ar- 
riving" at the state (1) as the result of collisions and second term describes solitons 
"leaving" this state. 

The general expression for the collision integral has the form: 

£ t = Y,^k{F t F k }. (5.7) 

k 

Here: 

C ik {F i F k } = J dl'd2'd2{W tk (l,2\l',2')F t (l')F k (2')- 

W ifc (l , ,2 , |l,2)F i (l)F fc (2)}. (5.8) 
Let us discuss ones more the detailed balance principle for usual particles in the form: 

W'(l / ,2 / |l,2) = W , (l,2|l , ,2 / ) (5.9) 

In this form detailed balance principle describes two particle collision (generalization 
to three, four, ... particle collision is well known). Formula (5.9) means that "arriving" 
number in state (l',2') from state (1,2) is equal to the "leaving" number from state 
(l',2') to state (1,2). If the total number of "arriving" particles to some fixed state 
(1,2) is equal to the total number of "leaving" particles from state (1,2) then new 
smoothened local balance principle can be formulated as: 
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J W ik (l, 2\1', 2')dl'd2' = J W ik (l', 2'|1, 2)dl'd2' (5.10) 
where i,k = s, b,ph. 

It is not difficult to show that the probabilities of collisions defined in section 3 
satisfy this condition. Thus the probabilities, W ik , of LE scattering processes in SG 
model satisfy the smoothened local balance principle smoothened local balance prin- 
ciple (5.10), while the dissipation parts of probabilities satisfy detailed balance 
principle (5.9). 

The distribution functions /, B, N in thermodynamic equilibrium must satisfy the 
following condition: 

C ik {F h F k } = 0, (5.11) 

in accordance with the smoothened local equilibrium principle for each i — k collision 
integral. It is necessary to emphasize that smoothened local balance principle (5.10) 
puts limitation on the probabilities of collisions, on the thermodynamic equilibrium 
condition (5.11) and on the distribution functions Fj. 

Having mentioned the general properties of collision integrals of LE in the SG 
equation, we proceed with the examination of each one. From formulas (4.21)- (4.22) 
it is easy to obtain the following expression for the collision integral for solitons in the 
low density case : 

Cs{h,F 2i } = -Sv s ^ + V s (v s )^ : (5.12) 

where renormalization of soliton velocity 8v and local coefficient of self-diffusion 
V s (v) are given by summing the partial contributions: 

5v s = ^25v si V s (v s ) =J2 T> si(v s )- (5.13) 

i i 

Here 

8v si = J \v \Ax si (vi s ,v 2 i)Fid2 (5.14) 

V *i = \j \v \[Ax si (v ls ,v 2i )} 2 F i d2 (5.15) 

The collision integrals for breathers can be analyzed by a similar way, but with 
a very important difference, connected with the conditions (4.20) and (4.21). The 
simplest case of breather ensemble is a ensemble with distribution function of the 
form: 

B(x, ip, v, uj 2 , t) = B(x, ip, v, t)8(u - u 2 ), (5.16) 

where Uo and density of particles rii satisfy the condition: cu m(v) » rii. 
It is possible to consider more general breather distribution functions, e.g. 

B = B(x, ip, v, t)b(u), 
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where b(u) has a sharp maximum near the point 00 = ujq. For simplicity only the case 
(5.16) will be considered. In this case all integrations under oo 2 are trivial. 
For breathers the collision integral has the following form: 

C h {B l , F 2i \ = -5v b — $u b —— + V b —— + Fh-^r + ^bjr^-- (5.17) 

ox dip ox 2 dip 2 oxoip 

The first two terms in this formula describe renormalization of breather velocity 5v b 
and its internal oscillation frequency 8uj b . The last three terms describe self-diffusion 
in (x, ip) space. As in the soliton case, the quantities 5v b , Sco b , T> b , T b and IC b are sums 
of partial contributions: 

i i i i i 

where: 

5v bi = J \v \Ax bi (v lb ,v 2i )F 2 id2 (5.19) 

5u bi = J \v \A(p bs (v lb ,v 2 ,i)F 2i d2 (5.20) 

V bi = \\ \v \[Ax bt {v lb ,v 2t )] 2 F 2i d2 (5.21) 

T ^ = \\ \vo\[^ bl (v lb ,v 2i )} 2 F 2l d2 (5.22) 

^ bi = \l \ v o\^ x bi(vib,V2i)Aip bi (v u ,v 2 i)F 2i d2. (5.23) 

6. KINETIC EQUATIONS AND ENTROPY PRODUCTION 

The Boltzmann type kinetic equations for LE with collision integrals constructed 
in previous subsection can be written following the standard way as: 

ft +<» + *»•<»»§ = ^w0 («> 

(6.2) 

diV (97V ON d 2 N d 2 N d 2 N 

-H^av))-^^))- = vAv)^™ A v)—^ A v) w 

(6.3) 
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Here terms from collision integrals describing velocity renormalization have been 
written on the left of eqs. (6.1)-(6.3). On the right of these equations there are only 
those terms describing dissipative processes leading to homogenization of distribution 
functions of LE. 

It is necessary to emphasize that collision integrals in kinetic equations (6.1) - (6.3) 
are equal to zero in the homogeneous case. Therefore the stationary solution of kinetic 
equations have following form: 

f = f(v), B = B(v,uj 2 ). (6.4) 

Here f(y), B(v,U2) are arbitrary functions of its arguments. 

In other words the kinetic equations (6.1)-(6.2) describe homogenization i.e. mixing 
of distribution function of LE up to homogeneous state in real space (for breathers - 
in (x, ip) space) and demonstrate that chaotization in momentum space cannot be 
realized. 

For chaotization processes in momentum space it is necessary to exceed the limits of 
exactly integrable model i.e. to take into account the terms destroying the integrability 
in the Hamiltonian of the system. 

Let us show now that homogenization of distribution function leads to entropy 
production in the SG localized excitations gas. The entropy of classical soliton gas and 
boson gases of breathers and phonons can be defined by standard way: 

S = ^2S k , k = s, b,ph (6.5) 

k 

S k = - f F k ln(F k /e)dl. (6.6) 
The entropy evolution in time is described by the formulae: 

dS k f dF k (l) , v , f . 

— = -J- Br ]nF k (l)dL. (6.7) 

Using kinetic equations (6.1)-(6.3), definitions of V, JC and T one can find that: 

dS k 



= J q k dl, (6.8) 



dt 

where the source q of the entropy production is: 

Qk = /EMl,2)|^{[Ax(l,2)] fel ^^ + [Av,(l,2)] fc ,^l} 2 rf2. (6.9) 

It is obviously that the expression is positive. It means that (6.9) proves the Boltz- 
mann entropy production theorem. We would like to emphasize that the entropy 
production is connected only with inhomogeneity in real space (x,</?).It is easy to see 
that in homogeneous case dF/dx = dF/d<p = and there no entropy production. 
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7.TRANSP0RT EQUATIONS 



In the present subsection the consequences from kinetic equations (6.1)-(6.3) will be 
analyzed. Let us emphasize that the homogenization of LE in real space means the 
homogenization of LE in temperature, concentration and macroscopic velocity spaces. 
The local macroscopic temperature T(x) is defined trough the local energy E(x) aver- 
aged over a distance do around x satisfying the inequality |AXj| << d . In other words 
collisions of LE lead to creation of diffusion, thermo conductivity and intrinsic friction 
processes. It is easy to note that the transport equations will have the form of local 
conservation laws for each type of LE separately. In Sine-Gordon system the numbers 
of solitons, breathers and phonons are conserved separately. Besides, momentum, en- 
ergy and angular velocity |^ (in the breather case) of each LE are conserved in each 
collision. The transport equation can be written in the following general form. 
For the solitons: 



d_ 
dt 



n s <a s >+ — [U: + Ur}=0. 



(7.1) 



For the breathers: 




(7.2) 



In formulae (7.1)-(7.2) following notation have been used. 
For solitons: 




(7.3) 




(7.4) 




(7.5) 



For breathers: 



nb < ab >= J a(x, v, if, u)B(x, v, ip, cu, t)dvdu 
Ul — / a(x,v, (f,u)[v + 5v]B(x,v, (p,u,t)dvdu 



(7.6) 



(7.7) 






(7.8) 
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W b r = J a(x,v,ip,uj)[uj + Suj]B(x,v,ip,uj,t)dvduj (7.9) 

d_ 

dtp 



d f 

W b m = — — — / a(x,v,ip,u)J r i ) B(x,v,(f,u,t)dvduj— 



d f 

— — J a(x,v,ip,u)lCbB(x,v,ip,Lj,t)dvdu}. (7-10) 

Substituting 1 , velocity and energy of LE for the variable a it is easy to obtain fol- 
lowing transport equations: continuity equations, hydrodynamics equations and equa- 
tions for local energy density. Let us write these equations in the explicit form. 

Continuity equations (a = 1) can be written as: 

Here the standard notation ji an d H have been used for Ui and Wi with a = 1 
respectively. 

Hydrodynamics equations can be derived from equations (7.1), (7.2) with a — v and 
(v,U2b) for solitons and breather respectively and have the following form: 

J^ sMs + JV; + Pf) = 0, (7.13) 

jn b u h + ^(Pl + PD + -^(Tll + EH = (7.14) 

^n b w b + ^(Ql + Qf) + -^(Rl + K) = 0. (7.15) 

Here u s and u b are hydrodynamic velocities of solitons and breathers correspond- 
ingly, w is hydrodynamic velocity in tp space, P[ and P™, % — s,b are pressures for 
solitons and breathers, the quantities HI and are pressures of breathers due to 
inhomogeneities in ip space, values Q and R are defined by formulas (7. 7)- (7. 10) with 
a = uj 2 b- 

The energy transport equations can be derived when a — E i for solitons and 
breathers. These equations can be written as: 

|n s T s + ^(f/; + f/7)=0, (7.16) 

^n b T b + JV[ + V?) + A W + W b m ) = 0. (7.17) 
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The quantities T{ means the average energies of corresponding LE. When Ui = Wi = 
the quantities Tj are the average energies of chaotic motion. Ui and Wi are energy 
density currents; one can conclude that only U™ and W m are different from when 
Ui = Wi = 0. 

Let us emphasize the important property of transport equations (7. 1 1)-(7. 17) . It is 
easy to see that for any homogeneous distribution functions: 

f = f(v,t), B = B(v,uj 2b ,t) (718) 

with constant temperatures, hydrodynamic velocities and chemical potentials /ij then 
all dissipative terms in these equations are equal to zero. In other words there no energy 
and momentum exchange between homogeneous gases of solitons and breathers. 

This special property of equations (7.14) - (7.16) is eliminated by taking into con- 
sideration terms in Hamiltonian which destroy the integrability of the model. 

As an example of explicit calculation of transport coefficients we shall obtain the 
expression for the self-diffusion coefficients assuming that the distribution function of 
solitons and breathers have following form: 

f = C s e ^bt, B = C b e - u ). (7.19) 

Let us discuss the concrete expressions for solitons and breathers diffusion coeffi- 
cients due to soliton-solitons, breather - breathers and solitons breathers collisions. For 
simplicity we will consider the nonrelativistic solitons and breathers only (v << 1). 
This case corresponds to small temperatures T « m. 

The diffusion current of solitons can be presented as: 

Js = + U b D sb] + -Q^ n s D sb (7-20) 

Using formulae for (Aip)^ and (Ax)ik from paragraph 3 , formulae (5.13), (5.18), 
(7.5), and (7.19) it is possible to calculate both D ss and D sb . We will present here the 
final results: 

D ss = (l/4)5 s 2 (T/ 7 rM s ) 1 / 2 {[ln( 7 M s /T)] 2 + C} (7.21) 
D sb = (l/2)5 2 (2T/7r^) s 1 6 /2 J sb , (7.22) 



where: 



I sb = (\n 1 -±^) 2 , if Tjix sb « 1 - ul 

1 — C^o 

I sb = [\n{2 1 u sb /T)f + C if 1 > T/fi sb > 1 - cu 2 . (7.23) 



Here: 
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C « 1.6, 7 « 1.8. 

There are two dissipative currents for breathers j^and if, which can be written as: 

diijf) $77 $7%} 

j* = (n b D bb + n s D bs )— + n b D bs —^ + (n b K bb + n s K bs ) 



r\ i ' "W" US r\ 1 V <J UU ' ' v s US I r\ 

OX OX 0(f 

$71/ u $71/ $71/ f) 

if = (n b F bb + n s F bs ) — + n b K bs -^- + (n b K bb + n s K bs )— (7.24) 

After routine calculations the coefficients D ik , F ik and K ik can be presented in 
following general form: 



D, 



- 5 1(JLYjd 



i 



D bs = -p^= (—) 2 J b D s , D ik = {D ik , F tk , K ik } (7.25) 

4V27T \Hbs. 



Here: 



M % M k . 

^ = MTM4' a = S ' 6 (7 ' 26) 
where M s and M b are defined by formulae (2.16). 

The expressions for Jf k are presented in tab lei. It is easy to see that the diffusion 
of breathers in real space ( x - space) is much faster then the relaxation on <p. 

ACKNOWLEDGMENTS 

Authors are grateful for Prof. R.N.Gurgi and Prof. L.N.Pastur for discussions. Part of 
this work was supported by INTAS grant Ref. No 94-3754. IVB and VGB are grateful 
for the hospitality of the Research Center of Crete/ FORTH. 



21 



References 



1. G. Zaslavskii and R.Z. Sagdeev, Introduction to Nonlinear Physics (Nauka, Moscow, 
1988). 

2. N.J. Zabusky and M.D. Kruskal, Phys. Rev. Lett. 15 , 240-243 (1965). 

3. C.S. Gardner, J.M. Greene, M.D. Kruskal and R. Miura, Phys. Rev. Lett. 19, 
1095-1097 (1967) 

4. V.E. Zakharov, S.V. Manakov, S.P. Novikov and L.P. Pitaevskii, Theory of Solitons 
(Moscow, Nauka, 1980) . 

5. L.D. Faddeev and L.A. Tahtadzhan, Hamiltonian methods in the theory of solitons 
(Nauka, Moscow, 1986). 

6. V.E. Zakharov, Zh.Eksp.Theor.Fiz. 60, 993 -1000 (1971). 

7. T.R. Koehier,A.R. Bishop, J. A. Krumhansl, J. A. and J.R. Schrieffer, Sol. State. Comm., 
17, 1515-1519 (1975) 

8. Y. Wada, J.R. Schrieffer, Phys.Rev. B 18, 3897-3912 (1978). 

9. N. Theodorakopoulos, Z. Phys. B 33, 385-390 (1979). 

10. K.Z. Fesser, Z. Phys.B 39, 47-52 (1979). 

11. K. Sasaki and K. Maki, Phys. Rev.B 35, 257 - 262 (1987) 

12. I.V. Baryakhtar, V.G. Baryakhtar and E.N. Economou, Phys.Lett. A207, 67-71 
(1995). 

13. V.G. Baryakhtar, I. V. Baryakhtar, B. A. Ivanov and A.L.Sukstansky, in: Proc. of the 
II Int. Symp. on Selected Topics in Statistical Mechanics, JINR, Dubna,417-426 
(1981). 

14. O.M.Braun and Yu.S. Kivshar, Phys.Rep. 306,Nl-2, 1-108 (1998). 



22 



Table!. Diffusion coefficients for breather 



ik 


T/mk < 1 - 
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r ln 1+^012 

1- 1— 


vrln ] +w " 
l-w n 
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bb 
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